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ABSTRACT 

We present results from a blind, spectroscopic search for redshift 5.7 Lya emission- 
line galaxies at Keck I. Using a band-limiting filter and custom slitmask, we cover the 
LRIS detector with low resolution spectra in the 8100 — 8250 A atmospheric window 
which contains no bright night sky emission lines. We find nine objects with line fluxes 
greater than our flux limit of 6 x 10 -18 ergs s _1 cm -2 in our ~ 5.1 square arcminute field. 
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We rule out a Lya identification for six of these based on the absence of the continuum 
break, expected at rest-1215 A for high-z galaxies, and/or the identification of additional 
emission- lines in our follow-up spectra. We find that extremely metal-poor, foreground 
emission-line galaxies are the most difficult type of interloper to recognize. For the three 
remaining emission-line objects, we identify a plausible counterpart for each object in 
a deep V-band image of the field suggesting that none of them has a continuum break 
in the i band. Our preliminary conclusion is that our field contains no z=5.7 Lya 
emitters brighter than 0.6L^ ya , where L* L = 3.26 x 10 42 ergs s _1 . Selecting a field 
with zero Lya emitters is marginally consistent with the no-evolution hypothesis - i.e. 
we expected to recover 2 to 3 Lya emitters assuming the Lya luminosity function at 
redshift 5.7 is the same as it is at redshift 3. Our null result rules out a brightening 
of L* Lya by more than a factor of 1.7 from redshift 3 to redshift 5.7, or, over the same 
redshift interval, an increase of more than a factor of 2.2 in the number density of 
Lya emitters. The paucity of z=5.7 Lya emitters raises the question of whether the 
Lya -selected population plays a significant role in maintaining the ionization of the 
intergalactic medium at z = 5.7. We find that if the escape fraction of Lya radiation 
is less than 0.4fL y c, where fi y c is the escape fraction of Lyman continuum photons, 
then the star formation rate in the Lya emitting population is high enough in the 
no-evolution model (our upper limit) to maintain the ionization of the IGM at z=5.7. 

Subject headings: galaxies: formation — galaxies: evolution — galaxies: high-redshift 
— techniques: spectroscopic 

1. Introduction 

The Epoch of Reionization likely marks a distinct change in the galaxy population in addition 
to a phase transition in the intergalactic medium. Sustained star formation is probably prevented 
by supernova feedback in the first, low mass galaxies. Completion of Reionization likely follows 
the formation of fairly massive galaxies that are immune to such violent feedback. The WMAP 
satellite has measured a large optical depth to electron scattering after cosmological recombination 
implying significant reionization at z ~ 17±5 (Kogut etal. 2003; Spergel etal. 2003). Yet the most 
recent re-heating must have occured more recently based on the thermal history of the intergalactic 
medium derived from the Lya forest (Wyithe & Loeb 2003). The absence of neutral regions in 
the intergalactic medium at z< 6 (Becker etal. 2001; Djorgovski etal. 2001) indicates Reionization 
was completed by this time. The population of galaxies responsible for Reionization, however, has 
yet to be recognized. 

Searches for z ~ 6 objects latch onto either the Lya emission line or the Lyman continuum 
break. Over cosmological distances, the HI opacity of the IGM blankets the continuum from 
912 A to the Lya line effectively shifting the continuum break up to 1216A (Madau 1995). For 



-3- 



galaxies at redshift 5 to 7, this break crosses the r and i bands, producing very red r — i galaxy 
color (or, as redshift increases, i — z color). Selection by very red color, i.e. r-band (or i-band) 
dropouts, provides many high-redshift candidates (e.g. Yan et al. 2003) but is not a sufficient 
selection criterion. Galactic L and T dwarf stars, as well as z ~ 1 galaxies with a 4000 A continuum 
break, have similar red colors in these bands. Broad-band selection followed by spectroscopic follow- 
up for Lya emission has provided some of the best statistics on high-z galaxies to date (Lehnert 
& Brewer 2003), but the spectroscopic follow-up required for confirmation is difficult. Indeed, to 
date, the highest redshift galaxies have been discovere fd in narrow-band imaging surveys for Lya 
emission (z=6.56, Hu etal. 2002; z=6.541, 6.578, Kodaira etal. 2003). Selection by Lya emission, 
however, is clearly not complete as only ~ 25% of starburst galaxies at z ~ 3 show strong Lya line 
emission (e.g. Steidel etal. 2000; Shapley etal. 2003). Since line selection and continuum-break 
selection have different shortcomings, progress will likely be made by using the two techniques in 
parallel. 

At z> 5, a method for detecting intrinsically faint galaxies - i.e. typical galaxies rather than 
rare objects - in large numbers is clearly needed to characterize their number density and clustering 
properties. Surveys for lensed objects along cluster caustics probe extremely deep but search tiny 
volumes (e.g. Ellis etal. 2001). In contrast, narrowband imaging covers large areas on the sky but 
does not go as deep as desired. The survey described in this paper explores a technique that may 
allow us to probe further down the luminosity function than narrowband imaging surveys while 
covering more area than lensing surveys. This multi-slit windows technique minimizes the sky noise 
under an emission line by dispersing the light but covers significant area on the sky by stacking 
multiple longslits on a single mask. We describe the strategy and our pilot observations in §2 of this 
paper. In §3, we present our catalog of emission-line objects and discuss the line identifications. 
We constrain the number density and evolution of high redshift Lya emitters in §4. Section 5 
summarizes our main results and the utility of the multi-slit windows technique. We adopt an 
VL m = 0.3, Qa = 0.7, h = 0.7 cosmological model throughout this paper. 



2. Emission-Line Search 

To search for Lya emission line objects at redshift 5.7, we have carried out a multislit windows 
search at Keck I with LRIS (Oke et al. 1995) over the atmospheric window at 8200 A. Dispersing 
the light with a grating maximizes the line vs sky contrast. A narrowband filter restricts the 
spectral bandpass to the airglow-free region illustrated in Figure 1. Because the spectra are only 
~ 100 A wide, we can cut multiple long-slits on a single slit-mask, thereby greatly increasing the sky 
area surveyed. Similar techniques have been described by Crampton and Lilly (1999) and Stockton 
(1999) in the atmopheric window at 9200AOur aim was to systematically probe the 8200 A window 
for Lya emission at z=5.7 going significantly deeper in intrinsic luminosity than the Crampton 
&; Lilly work and narrowband-imaging surveys. Two key features of our strategy were (1) near 
real-time spectroscopic follow-up of emission-line candidates, and (2) a quantitative analysis of our 
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selection function based on simulations of our detection efficiency. 

2.1. Experimental Technique 

The shape of the filter transmission profile is an important aspect of the experimental design. 
To optimally use the detector area, the filter transmission curve should have a top-hat profile 
with width less than or equal to that of the atmospheric bandpass. Figure 1 shows the filter 
transmission curve of NB8185 which Esther Hu kindly loaned to us. The transmission profile is 
centered at 8197 A and has a full width at half maximum intensity (FWHM) of 106 A and full width 
at 10% transmissivity of 145 A. The filter wings transmit bright sky lines, so we spaced the spectra 
165 A apart to avoid contaminating the targetted bandpass with sky lines from neighboring slits. 
The increased slit spacing, over that allowed by a perfect top-hat transmission profile, effectively 
reduces the survey area by a factor of 1.56. 

We designed our mask for the Keck I LRIS instrument and the NB8185 filter. Our selection of 
the 150 mm -1 grating (4.7 A pix -1 ) reflects a tradeoff between maximum sensitivity (somewhat 
higher spectral resolution) and survey volume (more slits at lower resolution). Rest-frame line 
widths of redshift 5.7 Lya emitters are 260 km s _1 (Bunker etal. 2003) to 340 km s _1 (Ajiki 
etal. 2002), so the ideal spectral reslution for our experiment is 7 to 10 A. We chose a slit width 
of 1'.'5 which projects to 4.7 pixels and yields a spectral resolution of about 22 A for a source filling 
the slit. At this dispersion, we fit thirty-three slitlets onto the mask at 165 A spacing. Figure 2 
shows the mask layout. Each of the longslits projects onto the full T.27 length of the Tektronix 
2048 x 2048 CCD. The combined slitlets subtend 5.09 square arcminutes on the sky. Stability bars 
prevent the thin aluminum mask from buckling; and small holes allow accurate positioning of the 
mask using foreground stars. The mask design does not cover the full LRIS field due to vignetting 
in the spectrograph, and the spectra cover just over 50% of the detector width in the dispersion 
direction. 

The effective search area for faint objects, however, is significantly lower than the area sub- 
tended by the slits. Partial occultation by the mask renders objects with fluxes near the limiting 
flux of the survey too faint to detect unless they land in the center of a slit. These slit losses 
are derived quantitatively in §2.3. The depth of our search along the sightline, set by the filter 
bandpass, also has a slight dependence on flux since we measure the product FT(X) where F is the 
line flux and T(A) is the filter transmission at the wavelength of the line. We include this effect in 
our modeling of the source counts. The sky noise (and therefore our sensitivity) was found to be 
fairly uniform across a bandpass equal to one filter FWHM. For a 1'.'5 (or l'.'O) effective slit width, 
our survey volume is 1100 (or 730) Mpc 3 . 

We choose a field centered at RA = 0:26:35.86, DEC = 17:11:35.15, J2000.0 This field covers 
the core of the rich cluster CL0024+16, redshift 0.397 (Czoske 2001). The field was chosen because 
it could be observed at a low airmass from Hawaii, the reddening is low Eb~v = 0.057 (Schlegel, 
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Finkbeiner, & Davis 1998), and deep broadband imaging had been obtained previously for the 
galaxy cluster. By selecting a cluster field, we also improve the chances of serendipitously finding 
a lensed object. 

2.2. Observations 

We obtained the search spectra 2002 November 6 with the Low-Resolution Imaging Spectrom- 
eter (Oke etal. 1995) on Keck I. A total of 21,600 s of integration time was obtained under clear 
skies. A series of exposures were made to allow us to track flexure, remove cosmic rays, and apply 
a 4" dither along the slits. Mask alignment was verified on alternating exposures by checking the 
flux of the alignment stars. The median seeing measured from the spatial profiles of the alignment 
stars was 0'.'8. 

Calibration followed standard procedures using IRAF 5 . The mean bias level of each frame 
was measured for each amplifier and subtracted from the detector section read by that amplifier. 
Pixel-to-pixel variations in gain were flattened by dividing each data frame by an appropriately 
normalized exposure of a quartz lamp. The images were aligned using integer shifts to compensate 
for instrument flexure and dithering. The 12 frames were averaged using cosmic ray rejection. The 
dispersion and spectral coverage were verified using Ar lines in an exposure of a calibration lamp. 
At such low resolution, it was not possible to directly measure the spectral resolution as all lines 
are blends. 

We derived the flux calibration from observations of spectrophotometric standards made with 
a l'.'5 wide longslit. Despite the presence of an atmospheric absorption band in the filter bandpass, 
the two stars we observed, Feige 34 and Hiltner 600 (Massey & Gronwall 1990), gave remarkably 
consistent values for the sensitivity (within a few percent). The systematic errors introduced by 
the unknown location of a source within the slit are much larger and are discussed later in this 
paper. During the portion of the night with no moon, the measured sky brightness in our bandpass 
was 3.7 x 10~ 18 ergs s _1 cm~ 2 A -1 . 

Cleaning the image with sky and continuum subtraction was necessary to make automated 
searching possible. Although the brightest atmospheric emission lines have been blocked out of our 
spectra, we subtracted local sky in the combined spectral image by using a running sky window that 
subtended a region from 1'.'06 to 3'.'20 in both directions along the spatial axis. The resulting S/N 
in the sky subtracted frame reaches the theoretical limit set by the shot noise from sky photons, 
so the amplitude of the random noise fluctuations varies with wavelength as the square root of the 



5 The Image Reduction and Analysis Facility is written and supported by the IRAF programming group at the 
National Optical Astronomy Observatories (NOAO) in Tucson, Arizona. NOAO is operated by the Association 
of Universities for Research in Astronomy (AURA), Inc. under cooperative agreement with the National Science 
Foundation. 
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intensity of the sky spectrum. We fit the continuum line by line for each panel of the sky-subtracted 
image and subtracted the fit. Emission-lines can be picked out of this cleaned image in a manner 
similar to object detection in direct images. 

2.3. Simulations 

We used simulations to gauge the sensitivity of our emission line search. We inserted artificial 
emission lines into our stacked and cosmic-ray cleaned - but yet not sky- or continuum-subtracted 

- data frames. We used artificial emission lines with Gaussian spatial and spectral profiles. The 
three combinations of parameters were: (FWHM SJ)atiai , FWHM A ) = (0.85", 18.4A), (1.07", 23A), 
and (1.28", 30 A). These choices were motivated by the expected velocity dispersion - rest-frame 
220 km s" 1 FWHM (observed line width of 6 A) and our spectral resolution - 22 A for an object 
filling the slit. Since more compact objects are easier to detect, these cases include the limit of 
an unresolved source but also explore the sensitivity to source profile. For each line shape, we ran 
simulations for a range of line brightnesses from 75 to 1000 digital units (which corresponds to 
2.91 x 10~ 18 to 3.89 x 10~ 17 ergs s" 1 cm" 2 for no slit loses). We avoided placing artificial emission 
lines on bright continuum objects, just as we avoided regions around such objects in our automated 
search. 

We automated the recovery of these artificial sources by tuning the parameters in the SEx- 
tractor software (Bertin & Arnouts 1996). We choose regions with at least 4 connected pixels with 
S/N exceeding 0.5 times the rms background fluctuations. The noise level was computed in a ring 
around a candidate, so it sampled both the sky and continuum noise. This approach proved more 
robust than the human eye which may be drawn to regions of high shot noise where bright sky 
lines intersect continuum sources. The sky aperture underestimates the noise in regions near the 
edges of the bandpass, where the sky rms was higher due to the presence of airglow lines, so we 
excluded these regions from the automated search. The fraction of artificial emission line objects 
recovered determines the completeness of our emission line search. For the Of! 8 FWHM, 18.4A 
FWHM sources, SExtractor finds nearly all lines brighter than ~160 DN, or equivalently, to an 
emission line strength of 6 x 10 -18 ergs s _1 cm -2 for a source in the center of a slit. Our results 
are summarized in Figure 3 which shows the completeness is higher for more compact lines. 

Because of slit losses, the position of an object within the slit affects its observed flux and 
hence impacts our ability to detect it. Consequently, the effective area of our survey depends on the 
intrinsic flux of the object: the faintest objects can only be detected if they fall near the slit centre 
while intrinsically brighter objects can be detected closer to the slit edge and hence over a larger 
area. Figure 4 shows how the effective area of our survey grows as the flux of the target increases. 
This plot includes the small, but measured, reduction in search area caused by bright foreground 
objects that fall on a slit. The fluxes of sources that are only detectable exactly in the slit center 

- where the survey area tends to zero - are 4.5 x 10~ 18 ergs s -1 cm -2 , 6.2 x 10~ 18 ergs s -1 cm -2 , 
and 8.0 x 10~ 18 ergs s -1 cm" 2 for our three fiducial object sizes. For objects that fall in the central 
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two-thirds of a slit (i.e. l'/O), the area searched is 3.3895 square arcminutes and the flux limit is 
6 x 10~ 18 ergs s _1 cm" 2 , 8 x 10" 18 ergs s" 1 cm" 2 , and 10 x 10~ 18 ergs s _1 cm -2 . In Section 4 
we will use the flux-dependent survey area function shown in Figure 4 to predict expected source 
counts. 

3. Object Identification 

The dividing line between faint sources and noise fluctuations becomes increasingly fuzzy as 
line flux decreases. Based on the analysis of our simulations (see §2.3), we believe most lines with 
160 counts or more are real detections. To ensure that we did not miss objects, we ran SExtractor 
(Bertin & Arnouts 1996) with the same detection criteria which were applied to the simulations 
in § 2.3 and identified emission lines objectively. Figure 5 shows the nine objects we found in the 
discovery spectra above this count limit. 

We extracted spectra for the nine emission-line objects, and Table 1 summarizes their measured 
line properties. Notice that several line sources - B5, B4b, A21, and C28 in in order of increasing 
flux - have line fluxes that if at z=5.7 would imply Lya luminosities similar to L* Lya , where 
L* Lya = 3.26 x 10 42 ergs s" 1 . We adopt this fiducial luminosity for the knee of the luminosity 
function based on the typical Lya luminosities measured at z=3 (Steidel etal. 2000) and at z= 
4.9 (Ouchi 2003). If L* Lya does not evolve, an L*^ ya emitter at z — 5.7 has an observed flux of 
1.09 x 10" 17 ergs s" 1 cm" 2 . Our faintest emission line source, B5, would have a line luminosity of 
0.77L* if it were at z = 5.7. For completeness, we remind the reader that we do not know where 
in the slit a particular source fell so the line fluxes are technically lower limits. 

Any viable high-redshift Lya emitter will have at least a factor of 5 to 10 continuum break 
just blueward of the Lya line due to attenuation from intergalactic gas (Madau 1995). Figure 5 
reveals relatively bright continuum emission to both sides of the line in B5, D21b, and D22. We 
reject these three objects as Lya candidates on the basis of the absence of a continuum break; and 
we tentatively identify D22 as a blend of Ha and [Nil] 6584 on the basis of the line shape. 

For the other galaxies, the follow-up observations described in the next section were required 
to determine their identity, but we draw attention to a few interesting attributes derived from the 
discovery spectra first. For example, extremely high equivalent widths are often more easily fitted 
by the Lya line than any other emission line, and we find C17 and D21 have extreme equivalent 
widths. Specifically, using the continuum blueward of the C17 line, we measure observed equivalent 
width of 490 A. No continuum is detected near D21, and we find the equivalent width must exceed 
500 A. In addition to their high equivalent widths, objects C17 and D21 are the only two candidates 
not spatially resolved in / . / 8 seeing. 
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3.1. Follow-Up of Emission-Line Candidates 

3.1.1. Follow- Up Spectroscopy 

To ensure accurate alignment of our follow-up mask on the desired targets, we milled the 
follow-up mask using a modified version of the mask design file from our search mask. Conse- 
quently, the relative positions of set-up star boxes and those slit segments that contained the five 
candidates of interest were preserved, allowing us to position our follow-up mask's slits exactly on 
the same place in the sky as those in our search mask. To allow us to use higher dispersion and 
broader spectral coverage for the follow-up observations without generating overlapping spectra, slit 
segments containing no Lya candidates were removed from the mill file. This procedure allowed 
us to be confident that our targets were indeed in the slits of our follow-up mask. 

We obtained follow-up spectra of emission-line objects A21, B4b, B5, Dll, and D21 on 2002 
November 9 with LRIS. We inserted the 831 mm -1 grating blazed at 8200A to improve the spectral 
resolution on the red side (0.93A pix -1 ). Blue light diverted by the D680 dichroic was dispersed 
by a 300 mm -1 grism blazed at 5000 A. We removed the narrowband filter to increase the spectral 
coverage from ~ A6800 to the atmospheric limit at ~ A3200 in the blue and A6800 to A8700 in the 
red. 6 

Five hours of integration time were obtained under thin cirrus clouds and / /8 seeing fwhm. 
We guided in the i-band to keep the ~ A8200 emission centered in the slits. Since there is no 
atmospheric dispersion compensator on LRIS, the targets are only in the slits on the blue side for a 
few hours. The measured width of the night sky lines was 4.7 pixels, or 4.4 A. The image processing 
steps were the same as those described in Section 2.2. Figure 6 shows the emission line sources 
A21, Dll, and D21 were redetected in our follow-up spectra. Object B4b was not re-detected. 

3.1.2. Imaging Follow- Up 

We mapped the positions of the emission line candidates onto a direct image of the sky taken 
with LRIS. Since we do not know the position of an object within the slit, the error perpendic- 
ular to the slit direction is ±0'/75. We computed only the lowest order term for the coordinate 
transformation along the slit, and the average error for a set of bright, test objects was (y/55. The 
coordinates on the LRIS image were tied to a deep V-band image of CL0024V (Czoske et al. 2001). 
None of our emission line candidates were near bright stars or galaxies. Object D21 lies along a set 
of arclets near the cluster although the spectral follow-up identified the line as Ha emission from 
a foreground galaxy. Although our spectroscopic analysis yielded 3 emission lines which could be 
Lya at z=5.7, we find a faint V-band object within 2'/0 of the position of each. The small size 
of these galaxies is consistent with either irregular galaxies at z ~ 0.6 or dwarf irregular galaxies 



Exact coverage varies by object according to the slitlit position. 
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at z ~ 0.24 The faintest counterpart, related to C28, consists of two small clumps. The brightest 
counterpart, related to B4b, appears to be a spiral galaxy near, but not on, the slit. Deep i or 
z imaging of these fields would determine whether there is another object closer to our predicted 
position. Follow-up spectroscopy of each V-band counterpart would also quickly determine if it is 
the emission-line source. 

3.2. Conclusions on Line Identifications 

We draw the following conclusions about the identity of the emission-line objects in Table 1. 

• A21: Several emission lines were detected in the follow-up spectrum. We identify the discovery 
line as [OIII] A5007; and the emitter is at z = 0.656. See Figure 6. 

• B4b: Emission-line source not detected on follow-up mask. 

• B5: Lack of a continuum break across the emission line rules out Lya identification. 

• C17: If the emitter is at z=5.7, the line luminosity is brighter (and more rare) than the 
objects expected in our survey volume. We tentatively identify a V band counterpart for this 
object which, if correct, indicates the object is not at high redshift. We cannot rule out a 
Lya identification at this time however. 

• C28: Marginal continuum detection to both sides of the line leads us to suspect this object 
is not Lya . Continuum detection is too weak to rule out a Lya identification however. 

• Dll: Several emission lines were detected in follow-up spectrum. We identify the discovery 
line as [OIII] A5007 at z = 0.662. See Figure 6. 

• D21: This object lies near the cluster caustic. No continuum or line emission, aside from the 
single line at A8189.65, was detected from D21 in the red follow-up spectrum (A7080 — A8970). 
The line profile has a Gaussian FWHM of 3.7 A which rules out an identification as [Oil] 
3726, 29 at z=1.2. However, a series of emission lines identified as [OIII] 4959, 5007, [OIII] 
4364, and the Balmer series are seen in the blue spectrum. The high equivalent width line 
originally detected turns out to be Ha emission at z=0.248. The 3-sigma upper limit on 
the [Nil] 6583 line, which is not detected, is < 1% of the strength of the Ha line. Our red 
spectrum covered [Arlll] 7135, Hel 6678, [SII] 6731, [Nil] 6583, [OI] 6300, and He 5578. None 
of metal lines are strong in this obviously very metal poor galaxy; and the He I lines fell on 
strong night sky lines. 

• D21b: Strong continuum emission, but lack of continuum break rules out Lya identification. 

• D22: Strong continuum detected, but lack of continuum break rules out Lya identification. 
Second line redward of brightest line presents correct offset for Ha and [Nil] 6583 at a redshift 
z ~ 0.24. 
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Objects B4b, C17, and C28 are our only objects for which a Lya identification is plausible. 
Within the positional uncertainties, the position of each of these emission-line sources matches the 
location of an object detected in a deep V-band image. In summary, we urge caution in the use of 
single-line redshifts for high-redshift galaxies. Many of the faint emission-line objects that we found 
turn out to be foreground galaxies upon more detailed inspection. Our sample contains at most 
three Lya emitters brighter than 7.8 x 10 -18 ergs s _1 cm~ 2 , and none of these three candidates 
appears likely to be a genuine z=5.7 galaxy 



Since our survey sensitivity and search volume are well understood from simulations, our 
observations constrain the properties of the z=5.7 population of Lya emitters. The clean detection 
of objects B5 and C28 confirm that a Lya emission line with L* luminosity would easily be detected, 
and we are confident that galaxies as faint as 0.6L* would be detected > 85% of the time. The 
paucity of Lya emitters found excludes a high density of bright galaxies with high Lya escape 
fractions. We quantify this statement in §4.1 and §4.2. In §4.3, we discuss whether bright, Lya 
selected galaxies can maintain the ionization of the intergalactic medium at z=5.7. 



To gain some insight into the properties of the high-redshift galaxy population, we parameterize 
the Lya emitting population by a Schecter luminosity function <f)(L/L*)d(L/L*) = 4>*(L/L*) a e~ L / L> d(L/L*), 
where the three parameters (ft*, L* , and a describe the number density, Lya luminosity of the 
bright-end turnover, and the faint-end slope. We expect our survey to find 



galaxies. The survey volume, V, and the survey completeness, £, depend on luminosity because 
the survey area and bandwidth decrease as the line flux approaches the limiting flux of the survey 
(see §2.3). We compute V{L/L*) from the cosmological volume increment per unit redshift, the 
accessible redshift interval (from Figure 1), and the effective survey area (from Figure 4). Figure 3 
illustrates the completeness function, £(L/L*). Guided by the properties of the Lya population at 
redshifts 3 and 4 (Steidel etal. 2000; Ouchi etal. 2003), we adopt the following luminosity function 

- 4>* = 0.0055 Mpc~ 3 , L* = 3.26 x 10 42 ergs s" 1 , and a = -1.2 (-1.6) - and integrate down to the 
limiting flux of the survey. For the no-evolution scenario, we predict our survey should return 2.5 
(2.4) Lya emitters per field on average. 

Using the above formalism, we test the hypothesis that the luminosity function does not evolve 
between redshift 3 and 5.7. Suppose we observe many fields, we would expect to recover the mean 

- i.e. 2 to 3 Lya emitters - on average. Assuming (for the moment) the galaxies are randomly 



4. 



Discussion 



4.1. Density of Lya Emitters 




(1) 
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distributed in space (i.e. Poisson statistics), the probability of obtaining (or 1) object is 8.2% 
(20%). The paucity of Lya emitters found in our field is therefore marginally consistent with 
no-evolution in the Lya emitting galaxy population between redshift 3 and 5.7. 

It is interesting to consider how strong lensing by the cluster CL0024+16 might affect the 
results of our survey. Galaxies seen through the cluster core will typically be magnified by a factor 
~ 2. The cluster core subtends ~ 18% of the area covered by our mask, so we would expect 1 in 5 
background galaxies to be lensed by a factor of two. We find no statistically significant clustering 
among the emission-line objects in our field. Only two of our emission line objects clearly have 
redshifts greater than the cluster redshift, and one of these, A21, is projected well beyond the 
cluster core. Given that we find no high redshift galaxies, accounting for the presence of a lensing 
cluster in our field would strengthen our claim that the emission line galaxy population does not 
undergo positive evolution from redshift 3 to redshift 5.7. Stronger magnification is only possible 
near the caustics where our mask subtends an insignificant area. 

Our results place strong constraints on the evolution of the Lya -emitting population. For 
any model of the luminosity function described by 4>*, L*, and a, we can calculate the number of 
Lya -emitters our survey would detect. Figure 7a shows the expected yield of Lya emitters per 
field over a grid spanning a factor of three evolution in L* and 4>*. For each of these models, we 
computed the probability of our field yielding zero Lya detections. The leftmost dashed line in 
Figure 7b illustrates the 99% Poisson confidence contour for our null result. We find that in the 
limit of no density evolution, the value of L* L appears unlikely to brighten by more than a factor 
1.8 between z=3 and z=5.7. Similarly, in the limit of no luminosity evolution, the number density 
of Lya emitters does not increase by more than a factor of 2.1 between z=3 and z=6. The number 
of predicted objects is not sensitive to the faint-end slope of the luminosity function because our 
observations reach only 0.6L*. 7 The other curves show the upper limits had we detected 1, 2, or 3 
Lya emitters. Even if a similar observation of another field yielded 2 Lya detections, we would 
conclude that between z=3 and z=6 L* brightens by less than a factor 3 and the comoving density 
of Lya emitters grows by less than a factor of 3.75. Our results rule out strong positive evolution 
in the Lya - emitting population from z=3 to z=6. 



4.2. Impact of Galaxy-Galaxy Clustering on Results 

The true variance in galaxy counts N about the mean fi will be larger than the Poisson variance, 
5 2 =< (N — n) 2 >= fi, by an amount reflecting the clustering strength, /U 2 £(r) (e.g. Peebles 1993). 
The correlation function, £(r), has not yet been measured for high-redshift Lya emitters, but we 
can illustrate how much clustering would be required to significantly change our conclusions. If 



7 Figuro 7a and Figure 7b use a faint-end slope of a — —1.2. For a steeper faint-end slope with a = —1.6, the 
allowed brightening of L* is slightly smaller (1.7 for zero Lya emitters detected) while the allowed increase in 4>* is 
slightly larger (2.2 for zero objects detected). 
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clustering reduces the number of Lya emitters in our field from 2.5 to 0, then Sp/p = 1 and 
£(r) = 1. The variance of the counts increases to 5 = \/p + P 2 C( r ) = 2.0, which is not much larger 
than the value of 8 = 1.6 expected for a random source distribution. Yet this estimate implies a 
higher degree of clustering than the Steidel etal. (1999) dp/p ~ 0.25 estimate for Lya emitters at 
z ~ 3. It seems unlikely, therefore, that clustering of Lya emitters could account for the paucity 
of Lya emitters found by our survey; and our conclusion - that L* Lya emitters cannot increase 
much in luminosity or increase in number density between z ~ 3 and z ~ 5.7 - holds even if the 
population is fairly strongly clustered. 

To place further constraints on the potential impact of clustering, it is interesting to ask 
whether our number density is consistent with the yield of high-redshift galaxies from other Lya 
surveys. For a non-evolving Lya population, we expected to recover only a couple Lya emitters in 
our pilot survey, and we expected them to have sub-L* luminosities. Narrowband imaging surveys, 
in contrast, cover larger areas but recover only objects brighter than roughly 2L*. The no-evolution 
model implies that a survey must, on average, cover 21 square arcminutes on the sky to return 
one 2.0L* Lya emitter from a 103 A wide bandpass at z=5.7. This area corresponds to a survey 
volume ~ 4500 Mpc 3 . The LALA group has three secure identifications of z = 5.7 Lya emitters 
which corresponds to 2.5 x 10 4 Mpc 3 per 2.5L* Lya emitter. For our adopted cosmology, the 
volume per emitter for the Hu et al. narrowband imaging is smaller at roug hly 6900 Mpc 3 . The 
Lya search at Subaru confirmed a bright z = 5.7 emitter (F = 4.1 x 10~ 17 ergs s _1 cm~ 2 ; Taniguchi 
etal. 2003) after searching about 3.4 x 10 4 Mpc 3 . Considering the uncertainties inherent in gross 
comparisons among these surveys, our limits on the number density of Lya emitters appear to 
be consistent with the other results. It remains plausible then that the Lya emitting population 
could maintain the ionization of the IGM at z = 5.7 and be related to the population responsible 
for cosmic reionization. 

Lehnert & Bremer (2003) selected galaxies between z=4.8 to z=5.8. using the broadband 
drop-out technique and spectroscopically confirmed two Lya emitters in the 8200 A atmospheric 
window. Their survey volume in the z=5.697 to z=5.710 window is only ~ 9.72 x 10 3 Mpc 3 , 
but both objects found are bright sources typical of those found by narrowband-selection (BDF 
1:10 F = 2.45 x 10~ 17 ; BDF 2:19 F = 2.51 x 10~ 17 ). 8 The implied number density of L> L* 
Lya emitters in the field is roughly 1 per 5000 Mpc 3 . Lehnert & Bremer (2003) conclude that 
their objects cannot maintain the ionization of the IGM and suggest that lower luminosity galaxies 
contribute to the ionization of the IGM. Since some of the emission line galaxies that we detected 
are too faint in the continuum to be selected by the Lehnert & Bremer survey, their result may 
be consistent with our conclusion. To date then, the different high redshift galaxy surveys seem 
to yield mutually plausible constraints on the high redshift galaxy population. Consequently, we 



They also find a much fainter source just redward of the window which has an observed equivalent width of 
40 A and would be difficult to recover in a narrow-band imaging survey. The object, BDF 1:19, has flux 3.09 x 
10~ 18 ergs s _1 cm~ 2 - just below the sensitivity limit of our survey. 
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cannot yet use the variance in Lya emitting galaxies among surveys to empirically estimate their 
clustering properties. 

4.3. Are There Enough High-Redshift Galaxies to Ionize the Intergalactic Medium? 

Recent observations of high redshift quasars pin the redshift of the final reionization epoch 
near z ~ 6 (Becker etal. 2001; Djorgovski etal. 2001), so the formation rate of stars and quasars 
must remain high enough following the end of Reionization to at least offset the recombination 
rate of the intergalactic gas. The identity of the objects that maintain the ionization of the in- 
tergalactic medium at z = 5.7 remains an open question. While quasars dominate the production 
of intergalactic Lyman continuum radiation by z = 0, at z ~ 3 more ionizing flux escapes from 
Lyman-break-selected galaxies than from quasars (Madau etal. 1999; Steidel etal. al 2001). The 
contribution from quasars at z ~ 5.7 appears to fall far short of that needed to maintain IGM ion- 
ization (Fan etal. 2001). It is interesting to examine whether the bright Lya emitting population 
is a plausible source of the required ionizing luminosity. 

The time from redshift z=5.7 to the reported completion of Reionization at z=6 (Becker etal. 
2001; Djorgovski etal. 2001) is only 152/i^y Myr. Hence the galaxy population responsible for 
Reionization may still be present at z = 5.7. A hydrogen atom in the IGM consumes several 
ionizing photons in this short period. The recombination time at z=5.7, averaged over a clumpy 
IGM, is t rec = 64 Myr (0.041/(J7fe/i7 ))C^ ) 1 , where C30 is the ionized hydrogen clumping factor 
in units of C =< n 2 HII > /n 2 Hn = 30 (Madau etal. 1999). To raise the filling factor, fun, of 
ionized regions in the IGM to unity, the number of ionizing photons escaping from galaxies over the 
recombination timescale must equal or exceed the baryon density. Following Madau etal. (1999), 
the critical production rate of ionizing photons is Nh = n#(0)/t rec (z) per unit comoving volume, 
which we can write as 

For an initial stellar mass function dM/dN oc M~ 2 - 35 (i.e. Salpeter IMF), continuous star formation 
produces ionizing photons at a rate of Njj = 10 53 ' 4 s" 1 per solar mass of 1.0 M to 100 M stars 
formed in 1 year (Leitherer etal. 1999). Including the mass of 0.1 to 1.0 M Q stars formed by this 
IMF, the minimum cosmic star formation rate required to keep the universe ionized at z=5.7 in 
the absence of any ionizing flux from quasars is 

A. = 0.0223 M yr^Mpc- 3 C3o (^) ' (^f ) ' f^, (3) 

where fi y c is the fraction of Lyman continuum photons escaping from galaxies. Theoretical sim- 
ulations predict that the escape fraction, fL y c, will become large when a starburst-driven wind 
blows out of a galactic disk (Fujita et al. 2003). In nearby starburst galaxies, measured values 
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of fiyC are less than 15% (Leitherer etal. 1995; Hurwitz etal. 1997). Some high-redshift galaxies 
have higher escape fractions (Steidel etal. 2001), but no clear consensus about escape fractions 
from high-redshift starbursts has been reached (Fernandez-Soto etal. 2003). Although the most 
appropriate value of fhyC for high-redshift galaxies clearly remains to be determined, values up to 
unity remain plausible for now. 

The Lya line is a notoriously poor measure of the galactic star formation rate. The large 
scattering cross section in the line produces a high probability that a Lya photon will be absorbed 
by a dust grain. The effective optical depth in the line can be dramatically reduced, however, by bulk 
motions in the ISM which Doppler shift the line photons away from the resonance. Observations 
of the Lya emission/absorption profile show wide variety among even starburst galaxies (Shaplcy 
etal. 2003), and models suggest the emission line flux depends as much on the evolutionary state 
of the starburst wind as on the absolute value of the star formation rate (SFR) (Tenorio-Tagle 
et al. 1999). For purposes of illustration, we simply parameterize the line luminosity in terms of the 
fraction of Lya photons that escape the galaxy fL ya - 9 For a Salpeter stellar IMF extending from 
0.1 Mq to 100 M a star formation rate of 1.0 M Q yr _1 yields an Ha emission line luminosity 
of 1.26 x 10 41 ergs s" 1 (Kennicutt 1998). Using the Case B ratio of the H recombination lines at 
low density (Ferland & Osterbrock 1985), we find Li ya = 8.32fL ya LH a . Substituting for the star 
formation rate in Equation 3, we derive the critical (comoving) luminosity density of Lya emission, 

W = 2.34 X 10- e rg Mpo-3 C X (l±i) 3 (20 , (4) 

required to maintain the ionization of the IGM at z=5.7. 

We may use Equation 4 to explore whether the population of z=5.7 Lya emitters can make a 
significant contribution to the required ionizing luminosity. For our fiducial model (0* = 0.0055 Mpc 
L* Lya = 3.26 x 10 42 , and a = —1.2), the Lya luminosity produced by galaxies brighter than 0.6L* is 
L = 0.0028L* Mpc" 3 . This population can only maintain the ionization of the IGM if fLya is small 
(relative to fiyc) - i- e - most of the emitted Lya is absorbed and therefore not detected. Using 
equation 3, our limits on the number density of galaxies brighter than 0.6L2 ya imply the escape 
fraction of Lya photons would have to be fi ya < 0. 39/i^c if these relatively bright galaxies were 
to maintain the ionization of the IGM. 

Since the optical depth in the Lya line is much higher than the optical depth at the Lyman 
limit, it is possible for a large fraction of the Lyman continuum photons to escape while most Lya 
photons are absorbed by dust after many scatterings. A Lya escape fraction fi ya ~ 0.4 also 
seems plausible when the properties of the known z ~ 6 galaxies are considered. For SSA22-HCM1 
(z=5.74), Hu et al. (1999) derive a SFR from the ultraviolet continuum that is twice as large as that 
derived from the Lya line. (Using the flux to SFR scale adopted in this paper, we get a SFR from 



Note that fhya is technically the escape fraction of Lya photons relative to the escape fraction of Ha photons, 
but we implicity assume all the Ha photons escape. 
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the continuum that is six times larger than that estimated from the line flux, so a value of f^ya in 
the range of 0.2 to 0.5 is at least plausible.) If the number of luminous Lya emitters at z=5.7 is 
similar to our inferred upper limit on their number density and the Lya escape fraction < 40%, 
then relatively bright Lya emitters with L ~ L* Lya can maintain the ionization of the IGM at 
z = 5.7. A significant deficit of L* Lya emitters relative to this no-evolution model would suggest 
that another population of galaxies, perhaps large numbers of dwarf galaxies, must be present to 
maintain the ionization of the IGM. That measurement should be possible soon. 

Continuum-break-selected surveys also find a paucity of galaxies at high redshift. In particular, 
Lehnert & Bremer (2003) conclude that the bulk of the ionizing flux that reionized the universe came 
from faint galaxies with Mab( 1700 A) > — 21. They find a lower density of z ~ 5.3 r-band dropouts 
than expected from no-evolution predictions tied to the lyman-break-selected galaxy population at 
z=3 and 4. The ultraviolet luminosity density from their z ~ 5.3 drop-outs is considerably less 
than that required to ionize their survey volume. Drop-out searches complement the Lya searches 
because the inherent biases are quite different. (The dropout technique yields some false-positive 
detections, but the Lya surveys may miss galaxies entirely.) Constraints from both techniques will 
probably be required to pin down the population producing the ionizing background at z ~ 6. 

5. Conclusions and Outlook 

We have shown that the multi-slit windows technique provides a very sensitive means of search- 
ing for emission-line galaxies. It probes significantly larger volumes than surveys of cluster caustics 
but reaches sub-L*objects which are largely inaccessible to narrow-band imaging surveys at high 
redshift. In our pilot survey, most of the emission-line sources were found to be foreground galaxies. 
It is possible that the Lya emitting population at z=5.7 is described by the same L* and 4>* as 
the Lya -emitting population at z ~ 3 to 4, and we recovered no Lya emitters simply because the 
richness of our field was poor compared to an average field. However, the paucity of Lya emitters 
in our field does strongly rule out significant positive evolution in the luminosity function. Neither 
the characteristic luminosity nor number density of Lya emitters can increase much between z=3 
and z=5.7 - a period long enough, l.lh^ \ Gyr, to allow several cycles of bursting activity and 
fading (Sawicki & Yee 1998; Shapley etal. 2001). 

Since the recombination time of the intergalactic gas is shorter at higher redshift, the paucity 
of Lya emitters at z ~ 5.7 raises the question of whether the Lya -selected population can keep 
the IGM ionized at z = 5.7. In the no-evolution scenario, in order for the inferred star formation 
rate to be high enough to maintain the ionization of the intergalactic medium, the fraction of Lya 
emission escaping must be less than 40%. This average Lya escape fraction does not appear to be 
unusual. Indeed, only 25% of starburst galaxies at z=3 have Lya in emission (Shapley etal. 2003). 
Our constraint is strictly an upper limit, however. Additional observations will either recover 
Lya emitters in number or push this upper limit low enough to require another population of 
galaxies (e.g. perhaps dwarf galaxies) to maintain ionization. The existence of another population 
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of dwarf starbursts is suggested by the detection of a z=5.7 Lya emitter with de-lensed flux of 
only 2 x 10 -18 ergs s _1 cm -2 in a search volume an order of magnitude smaller than that of our 
survey (Ellis etal. 2001). 

Future surveys with the multislit windows technique are expected to produce about a dozen 
high-redshift galaxies per pointing. The gains will come from a combination of improved sensitivity, 
the larger field of view provided by new cameras, and multi-night observing campaigns. Dispersing 
the spectra to a rest-frame velocity dispersion of 200 km s _1 , which is about ~ 5.5 A, will reduce 
the sky noise by another factor of 3 to 4. The throughput can probably be doubled by using a 
more efficient dispersing element. It should be possible, therefore, to reach fluxes a factor of 2.5 
to 3 times fainter than our limit of 6.2 x 10 -18 ergs s _1 cm -2 in a single night of observing time. 
Although the cost of higher dispersion would normally be reduced area, the new generation of large 
format CCD arrays more than makes up for the decrease in survey area. 

Atmospheric windows at longer wavelength allow the technique to be extended to higher red- 
shifts. However, even in the next window at 9200 A (z = 6.5), the exposure time required to reach 
a given luminosity is nearly a factor of two higher than that required at z=5.7 (in the 8200 A 
window). We would expect a large multislit windows survey in the 8200 A window to yield a secure 
measurement of the number density of dwarf (i.e. sub L*) Lya emitters at z=5.7. 

These observations would not have been possible without the narrowband filter kindly lent by 
Ester Hu. We thank David Crampton for his encouragement and useful discussions. CLM and 
MJS recognize and acknowledge the very significant cultural role and reverence that the summit of 
Mauna Kea has always had within the indigenous Hawaiian community. We are most fortunate to 
have the opportunity to conduct observations from this mountain. This work was funded in part by 
grants from the David and Lucile Packard Foundation (CLM) and the Alfred P. Sloan Foundation 
(CLM). 
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Table 1: Properties of Emission-Line Sources 


Object ID 


FWHM 


FWHM 


Line Flux 


Line Flux 


EW 


Line ID 


(MS-0hl7) 


(") 


(A) 


(Counts) 


(10~ 17 ergs s" 1 cm -2 ) 


(A) 


(redshift) 


A21 


1.47 


18.0 


340 ± 50 


1.30 ±0.2 


28 


[OIII], z=0.656 


B4b 


1.48 


15.1 


290 ± 29 


1.1 ±0.1 


> 80 


Unlikey Lya Line 


B5 


1.40 


10.0 


220 ± 50 


0.84 ±0.2 


9 


Not Lya 


C17 


1.10 


17.2 


1200 ± 44 


4.80 ±0.2 


490 


Unlikely Lya Line 


C28 


1.29 


19.5 


430 ± 56 


1.7 ±0.1 


80 


Unlikely Lya Line 


Dll 


1.08 


16.5 


640 ± 38 


2.5 ±0.2 


47 


[OIII], z=0.662 


D21 


0.86 


17.7 


6280 ± 66 


24.4 ±0.3 


> 500 


Ha at z=0.248 


D21b 


1.38 


20.8 


700 ± 90 


2.7 ±0.4 


32 


Not Lya 


D22 


1.16 


20.4 


2760 ± 110 


10.7 ±0.4 


40 


Likely Ha ±[NII] 



Table Notes - 

(1) Object name on Martin & Sawicki mask for field at Oh 17°. 

(2) Extent of emission line along the slit. 

(3) Observed line width. All lines are unresolved. 

(4) Counts in the emission line. 

(5) Line flux derived assuming the average filter transmission over the bandpass. 

(6) Emission-line equivalent width. 
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Fig. 1. — Relative intensity of the Mauna Kea night sky (courtesy of Tom Slanger). The transmis- 
sion of filter NB8185 (kindly loaned by Esther Hu) illustrates the atmospheric window probed by 
our survey. 
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Fig. 2. — Slitmask layout for the search mask. The mask consisted of 33 longslits interrupted 
by three bars needed for structural rigidity. Small squares are holes for alignment stars used for 
positioning the mask on the sky. Dispersion direction is left to right. 
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Fig. 3. — Survey sensitivity. The fraction of simulated objects recovered by the automated detection 
routine for the three combinations of object size parameters. We define our completeness limit at 
160 counts, which, in the absence of slit losses, corresponds to 6xlCT 18 erg s _1 cm~ 2 . 
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Fig. 4. — Our survey area as a function of emission line flux. The horizontal line denotes the 
geometric angle subtended by the slitlets. Slit losses reduce the transmitted flux, and the faintest 
emission lines are only detected when the objects land in the center of the slit. 



Fig. 5. — Emission-line candidates found by the automated search. The arrows denote 5" and 
100 A in the spatial and spectral directions respectively. Wavelength increases from left to right. 
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Fig. 6. — Follow-up spectra for objects A21, Dll, and D21. None of the confirmed line identifica- 
tions are Lya . 
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Fig. 7. — (a) Predicted survey yield for various evolution models of the Lya luminosity function. 
The axes indicate the amount of evolution in the number density (vertical axis) or luminosity density 
(horizontal axis) relative to our baseline luminosity function ( L* Lya (z = 3) = 3.26 x 10 42 ergs s -1 , 
4>* Lya {z = 3) = 0.0055 Mpc -3 ). The faint end slope of a=-1.2 was assumed for the luminosity 
function, but the results do not depend significantly on the assumed a. For a Lya luminosity 
function that does not evolve between z~3 and 2=5.7, we expect to recover ~2.5 z=5.7 Lya 
emitters on average, (b) Constraints on the evolution of the Lya luminosity function from z~3 to 
z=5.7. The leftmost curve — labelled "0" — shows the upper limit (Poisson 99% confidence) on 
the amount of luminosity and density evolution in the Lya luminosity function allowed by the null 
result of our survey. The area of parameter space to the right of the curve is ruled out by our lack 
of Lya detections. The other curves show the upper limits had we detected 1, 2, or 3 Lya emitters. 
Evidently between redshift 3 and 6, L* L does not brighten by more than a factor of 1.8 and ^>2i/a 
does not increase by more than a factor of 2.1. 
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